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Abstract 

In semiconductor heterostructure, the Coulomb interaction is responsible for the electric current 
drag between two 2-d electron gases across an electron impenetrable insulator. For two metallic 
layers separated by a ferromagnetic insulator (FI) layer, the electric current drag can be mediated 
by a nonequilibrium magnon current of the FI. We determine the drag current by using the semi- 
classical Boltzmann approach with proper boundary conditions of electrons and magnons at the 
metal-FI interface. 



PACS numbers: 72.25-b, 73.30.Ds 




FIG. 1: Schematics of the NM/FI/NM trilayer structure. 

The conventional Coulomb drag effect [TH3] occurs in two-dimensional electron gases sepa- 
rated by an insulator barrier. When one of the electron gas carries a current, the momentum 
transfer due to Coulomb interaction leads to a small current in the other electron gas. Re- 
cently, this current drag phenomenon has been discovered in a different system with entirely 
different physical mechanisms [4J: when an electric current is injected into a Pt bar de- 
posited on a magnetic insulator Yttrium-Iron-Garnet (YIG) film, it is found that a small 
electric voltage is induced in the other Pt bar, which is also deposited on the same YIG 
film but is located several millimeters away from the current carrying Pt bar. The authors 
(4j attributed their finding to the combined effects of spin transfer torque (STT) [fj E] and 
spin pumping [8] : the spin Hall [9] current generated by the electric current in one Pt 
layer (Pt is a known material with a large spin Hall angle) is absorbed by the FI and for a 
sufficiently large STT, the magnetic moment of the FI begins precessing. The precessing FI 
pumps out a spin current to the other Pt layer. Finally, the resulting spin current converts 
to an electron current due to the inverse spin Hall effect | [T0l - IT5] . 

In this Letter, we propose a different geometry in which the magnon current flows normal 
to the plane of the layers throughout the structure. We show that the electron spin current 
in the metallic layers induces a nonequilibrium magnon current in the FI layer. By using 
semiclassical Boltzmann approach for electrons and magnons, we are able to self-consistent ly 
determine these currents and thereby obtain the drag current for given geometrical and 
material parameters. The resulting drag current is several orders of magnitude larger than 
that in the nonlocal geometry in Ref. [4] 
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To be more specific, we consider a simple trilayer structure, shown in Fig. 1 schematically, 
where a ferromagnetic insulating (FI) layer is sandwiched by two heavy metal films (NM1 
and NM2) such as Pt and Ta. A charge current parallel to the plane of the layers is injected in 
the layer NM1. To determine the drag current in the layer NM2, we first establish transport 
equations for each layers and then find proper boundary conditions to solve the transport 
coefficients. 

Electron current and spin accumulation in Metallic layers. For the NM layers, a spin 
dependent Ohm's law has been well established and may be written in the following form 

j = |E+| (exct-ctxE) (1) 

where the spinor current density j and the electric field E are 2x2 vector matrices in spin 
space, a is a Paul vector matix, and c and are the electric conductivity and spin Hall 
conductivity respectively. The second term is the spin Hall current whose anti-symmetric 
form is essential for j to be an Hermitian in spin space (also noted that Excr^— <rxE 
due to non-communitivity of the Pauli matrices). The electrical field is related to the spinor 
chemical potential jl via E = — (l/e)V/i where e(< 0) denotes the electron charge. While 
it is possible to work with an arbitrary choice of the spin quantization, we proceed below to 
a special case where the magnetic moment of the FI is oriented in the z-direction and the 
electric current flows in the y direction. If we choose the spin quantization axis parallel to 
the z-axis, one can simply work on the two-component (spin-up and spin down) form of the 
Ohm's law; i.e., 

^) = ^)-^) ( 2 ) 

and 

where a = ±1 represent spin up an down. To determine the spin dependent electric field, 
we recall the spin diffusion equation [T7j 

^\x)-A*)\ = ^ {X) ~/ {X) (4) 

and its solution 

/J(x) - n l (x) = A ie - x/x °f + B ie x/X °f (5) 
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where X s f is the electron spin diffusion length, and the constants A4 and Bi (i = 1 for 
the NM1 layer and i = 2 for NM2) are determined by the boundary conditions. Although 
these equations apply to both NM1 and NM2 layers, they have different constraints set by 
experimental measurement. For the NM1 layer, we take Ey{x) = E ext where E ext is the 
applied electric field in the NM1 layer, while J NM2 dxj y (x) = in the open circuit of the 
NM2 layer. 

Magnon current and magnon accumulation in the FI layer. For the FI layer, we start with 
a general magnon Boltzmann equation in the presence of spatially dependent temperature 
T(x) and magnetic field H(x), 

dN m dN m dT dN m dH dN m (dN m \ 

+ V *^FTx + Vx ~dH ' ^ + = " {-W) scatt (6) 

where N m (x, q,T(x), H(x)) is the magnon distribution. The first term describes magnon 
diffusion. The second and third terms are responsible for the magnon transport in the 
presence of temperature and magnetic field gradients, which have been recently studied in 
the content of spincalorics fT8H2T] . The last term on the left side of Eq. (6) is associated 
with acceleration of magnons by external forces such as a confining potential at boundary 
[22]. The scattering term on the right side of the Eq. (6) may be modeled by the relaxation 
time approximation 

/ ON \ N - N N — N° 

\ 9t ) scatt. T m T th 

where N m (x) = J rfqA^ m (x,q)/ J rfq is the momentum averaged magnon distribution while 
A^(x, q) = eeq/fcg 1 T(x)_ 1 is the local equilibrium magnon distribution, where e q = Dq 2 + A g is 
the magnon dispersion, D is the spin wave stiffness, is the spin wave gap, and 
the x component of the magnon velocity. The first relaxation term describes those processes 
which conserve the number of magnons. For example, magnon scattering by a paramagnetic 
impurity has the form of V^a+a^] i.e., the impurity or surface roughness [23j [24] scatters 
the magnon q' to the magnon q. As long as we neglect the wave number dependence of 
the scattering matrix l/ qq /, this process can be modeled by the first term of Eq. (7). The 
second term of Eq. (7) does not conserve the number of magnons. The magnon absorption 
and emission relax the nonequilibrium magnons to equilibrium ones, e.g., magnon-phonon 
interaction [1]. 

For the present system, we consider uniform temperature and magnetic field, and there 
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is no external force on magnons. Then, Eq. (6) and (7) reduce to 

dN m (x, q) N m (x, q) - N m (x) N m (x, q) - (q) 

Vx o • l°J 

ox r m r th 

We may proceed to solve N m by the same way as for the electron distribution in mag- 
netic multilayers [T7j . Particularly, one may expand the nonequilibrium distribution by the 
Legendre polynomials, 



7V m (*,q) = A£(q)-' dN ° m{q) 



(9) 

n=l 1 

where /i m (x) is the n = component of the nonequilibrium distribution and 9 is the angle 
between q and x axis. By placing the above equation into Eq. (8) and by utilizing the 
orthogonality property of the Legendre polynomials, one can arrive at a series of algebraic 
equations for the coefficients g <<n \x). In the supplemental material, we show the solutions in 
some limiting cases. Once the distribution functions iV m (x,q) are determined, we can find 
the magnon accumulation and magnon current via 

jm(x) = J dqv x N m (q,x); 5n m (x) = J dq [N m (q,x) - A^(q)] (10) 

where \ib is the Bohr magneton. Note that a magnon carries spin moment —^fi{— —2jjL B ) 
where 7 is the gyromagnetic ratio. 

We may further simplify the solution of the non-equilibrium magnon distribution by 
discarding high orders (n > 2) of the polynomials. Consequently, we find a local relation 
between magnon accumulation and magnon current, 

^i m (*) = 2, B *==M (11 ) 
dx r th 

and 

jm(x) = l J^ 5n rn(x) (12) 

where I n are integration constants I n = / d 3 qv n dN ™^ . We point out that this local 
current expression is valid in the limit r t h 3> T m which is a good approximation for ferro- 
magnets [1] (see supplemental material). By combining Eqs. (11) and (12), we obtain the 
diffusion equation for nonequilibrium magnons, 

f- 2 Sn m (x) - = (13) 



where the magnon diffusion length is defined as l m = y ^^th^m- At room temperature 
(T = 300X), for YIG with ~ 10~ 5 eV, r m ~ 10~ 7 s and r th ~ 10 _6 s 5 l m is estimated at 
0.05cm, consistent with the measurement [26]. Equation (13) has the general solution, 

Sn m (x) = A F e~ x/lm + B F e x/lm (14) 

and thus the magnon current density reads 

jm ( x ) = (-A Fe -^m + BFe */im\ (15 ) 

Boundary conditions. The outer-boundary conditions at x = — L\ and x = d + L 2 , 
where Li, d and L 2 represent the thicknesses of the layers of NM1, FI and NM2, are 
jl^(-Li) = jl^(d + L 2 ) = 0. The boundary conditions at the metal-FI interfaces depend 
on the interaction between electrons and magnons. Here we assume a s-d type interaction 
— Jsd^ ' where a is the itinerant electron spin of the metal layer and is the local spin of 
the FI layer at the interfaces. The interaction conserves total angular momentum and thus 
the first boundary condition is the continuity of total spin current at the interfaces; i.e., 

(-/We) K(0") - ji(0-)] =Jm(0 + ) 

and 

jm(d-) = (-/We) [Jl(d + ) - Ji(d + )} (16) 

The total angular momentum current conservation simply states that electron spin current 
in the metals must be converted into magnon current in the FI layer at the interfaces. If 
the interfaces have magnetic roughness, the spin-flip scattering by magnetic impurities can 
transfer spin angular momentum to lattice via spin-orbit coupling. In this case, the outgoing 
spin current would be reduced [27] . 

The other boundary conditions at the interfaces should relate the electron spin accu- 
mulation to the non-equilibrium magnon density. Within the s-d model, one can treat the 
electron spin density as an effective magnetic field on the interface spin of the FI layer; i.e., 
H e ff = J s dSm z and thus 

/i t (0-)-^(0-) = ^n m (0 + ) 

and 

e8n m {dr) = ^{d + ) - fi\d + ) (17) 
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where 

e = HnDf/2 (18) 

J sd D(e F )a 3 -VkFTLt h (e-^BT) ^ > 

D (s F ) is the electron density of state at Fermi level, a is the lattice constant of the NM 
layer and Li s (z) = Y^k=i fi" ^ s the polylogarithm. The detailed derivation of e is arranged 
in the supplemental material. We note that Takahashi et al. have also proposed a boundary 
condition at the interface [28j which relates magnon spin current to spin accumulation; i.e., 
j m oc fjft — /A However, such boundary condition is unable to self-consistent ly determine 
the magnon current. On the contrary, the boundary conditions we have derived are able 
to uniquely determine the spin and magnon currents throughout the structure. A rough 
order of magnitude estimation of e can be readily obtained by using the following plausible 
parameters appropriate for Pt/YIG/Pt structure: D(e F ) ~ 3n e /2e j p, n e = 5 x 10 22 cm -3 , 
J s d 1 eV, ep = 5 eV, a = 4 i, and D = 6 meV • nm 2 , and thus e ~ 0.2 meV • nm 3 . 

With the above boundary conditions, the constants in Eq. (5) and (14) can be readily 
determined. If one uses an Ampere meter [29] to measure the average (measured) induced 
electric current density ffi = (l/L 2 ) J dxjy 2 \x) in NM2 layer, we find the ratio between the 
induced current and the injected current magnitude r] (= j^p /jy 1 ^ is, 



°h \ 2 ( A s / 



sech (ft) 


cosh 


- \ 




b~ l + 6tanh ^) 


sinh {€) + 


1 + tanh 


(ft). 


cosh (£) 



(19) 



where b = c£T t h/(4:e 2 l m \ s f), and L\ 3> A 5 / is assumed for simplicity. The first prefactor 
(ch/c) 2 origins from the two successive conversions between electric current and spin Hall 
current in NM1 and NM2 due to the spin Hall and inverse spin Hall effect respectively. The 
second prefactor A s //L 2 indicates that the maximum range of the current density in NM2 

(2) 

is A 5 /; i.e., if the thickness of NM2 exceeds A s /, the average current density jy would be 
inversely proportional to L 2 . Interestingly, when L 2 is much smaller than A s /, the induced 
electric current is also small; this is because the spin current at the surface x = d + L 2 is 
set to zero and thus the self-consistent calculation demands a small current throughout the 
NM2 layer. In Fig. 2, we show rj as a function of the thickness of the metal layer (NM2) for 
Pt/YIG/Pt trilayers with several different YIG layer thicknesses. We choose the material 
parameters as follows: Pt layer conductivity cp t ~ 0.1 (fid • cm) -1 , spin diffusion length 
X s f = 7 nm and the spin Hall angle c^jc = 0.05 [30]; magnon diffusion length l m = 0.05 cm 
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FIG. 2: The ratio of the average induced current density and the injected current density as a 
function of the NM2(Pt) layer thickness for three different thicknesses of the FI (YIG) layer. 

and magnon relaxation time r t h = 10 -6 s. We see that r] decreases as the thickness of the 
YIG layer increases due to the decay of magnon diffusion current. Also for fixed YIG layer 
thickness, rj reaches its maximum around L 2 = \ s f. The peak value of rj is of the order of 
10 -4 . If the injected current density is 10 6 A/ cm 2 , the induced voltage of a Pt bar with its 
length w = 1 cm would be = wj^p jc ~ 1 mV. If one replaces Pt with Ta which has 
larger spin Hall angle of 0.1 [31j, can be further increased by a factor of 4, which is 
rather significant and easily detectable experimentally. 

Finally we comment on the relation of our calculation with the experimental measurement 
|4j. In their experiments, when the first Pt layer injects a spin current to the FI layer, the 
magnons propagate in the plane of the FI layer in order to reach the second Pt. While 
there is a similar non-equilibrium magnon density buildup near the second Pt layer, the 
direction of the magnon current and the gradient of the magnon density are in the plane of 
the layer. In another word, there is neither magnon current nor magnon density gradient in 
the direction perpendicular to the layer such that the second Pt layer is unable to receive 
any spin angular momentum from the FI layer. Thus, we conclude that the nonlocal setup 
in the experiment [4J is not relevant to our theory. In the conventional nonlocal metallic 
spin valve, however, one does observe a voltage change of the entire detection bar due to 
the spin accumulation (not the spin current or gradient of the spin accumulation) in the 
channel. In the present case, we derive the induced current in the second Pt bar which is 
related to the spin current (or magnon density gradient) in the direction perpendicular to 
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the layer. Furthermore, the observed current in the experiment [I] has been attributed to the 
STT and spin pumping, which is several orders of magnitude smaller than what we predict 
in our geometry. 

This work is supported by NSF. 
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SUPPLEMENTAL MATERIALS 



A. Magnon diffusion equation in ferromagnetic insulator 



In this section, we derive the magnon diffusion equation, Eq. (13), from the magnon 
Boltzmann equation, Eq. (8). By placing Eq. (9) into Eq. (8), we have 



q dx de q \T m T th J 

<9A^(q) d\x m (x) dN^(q) \x m (x) \i m (x) 



oe 



Vr 



dx 



dN°Jq) dNl 



(Al) 



where = Jdq^/Jdq, N° m = J dqN° m (q)/ / dq, and g(x,q) = ta {n) {x)P n (cos 9) 

n=l 

is the Legendre polynomial expansion of the nonequilibrium magnon distribution when a 
rotational symmetry is assumed. 

The non-equilibrium magnon density and magnon current, as defined in Eq. (10), are 
simply the zeroth and first components of the Legendre polynomials, 



5n m (x) = AnI^^x) 



(A2) 



and 



where we have defined 



L = 



d 3 qv 



n dK(q) 



(A3) 



(A4) 



(2tt)3 J * de q 
v — (m/q) ' ^ e m(q)/^ is the magnitude of magnon velocity. 

The relations among can be readily obtained by multiplying Eq. (Al) by v™ and 
integrating over q, 

d ( nr.} 

(A5) 



-i-Jm{x) = 2fi B , [n = 0) 

dx T t h 



(A6) 



and 



L n+1 



n 



dg( n -V n + 1 dg( n+ V 



+ 



2n — 1 dx 2n + 3 dx 



+ /„(— + — ^ (n) =0, (n>2) (A7) 
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where we have used the following orthogonal relations for the Legendre polynomials 

1 

J duP n , (u) P n (u) = ^-^£ njn , (A8) 
-l 

and 

1 

/2 (n + 1) 2n 
du Pl (U) P n , (U) P n («) = (2n + 1)(2n + 3 /n> + l + (2n+1)(2n _ 1) ^,n-l (A9) 

-1 

At this point, the problem of solving the magnon distribution is equivalent to determining 
the infinite function series g^ n \x). To simplify the problem, we assume r t h 3> T m , which 
is generally valid for most insulating ferromagnetic materials [lj. In this case, the ratio of 
g(n+i) jg(n) j g ^out ^T m /T t h> To see it, we first neglect in Eq. (A6) and substitute the 
resulting expression for into Eq. (A5). One immediately sees that the spatial derivatives 
of both ji m and j m scale as l/y/r m T t h. Putting this scaling back to Eq. (A6) and (A7), we 
confirm that neglecting g^ (n > 2) is justified as long as r t h 3> r m . Combining Eqs. (A2), 
(A3), (A5) and (A6), we arrive at the magnon diffusion equation in the main text, Eq. (13). 

We emphasize that the validity of the magnon diffusion equation rests on the condition 
that the magnon conserving scattering is much stronger than the magnon non-conserving 
scattering, i.e., r t h 3> r m ; this condition is similar to the spin diffusion equation for electrons 
where the momentum scattering (spin-conserving) is stronger than the spin-flip scattering 
[2]. Thus, the magnon diffusion equation can be used for the length scale larger than the 
spin-conserving scattering length (e.g., it is the mean free path in the electron case) which 
is considered much smaller than the magnon (spin) diffusion length. 

B. Boundary conditions at metal-FI interfaces 



In this section, we derive the boundary conditions Eqs. (16) and (17) by using the s — d 
exchange coupling 

h = -j sd J2 1 d " rs w • s ( R *) 6 ( r - R ) > ( B1 ) 

i ^ 

where S(R$) are the localized spins of FI at the interface in contact with the conduction 
electron spins s(r) of the metal. The above s-d interaction contains a spin-flip scattering 
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processes of J s d(c^c^- C ^a (l + h.c.) and a non-spin-flip process J s dc£ + , CT Ck+q<7aqa q '|3j. For the 
spin-flip process, the electron spin loss (gain) leads to the magnon generation (annihilation), 
but the total angular momentum is conserved; this leads to the first interface boundary 
condition of Eq. (16). For the non-spin-flip processes, we use the simple mean field approach, 
i.e., the conduction electron spin density is treated as an effective magnetic field acting on 
the interficial spins of the FI layer; i.e., 

SH eff = — (// - (B2) 

where D{e F ) is the electronic density of states and v is the volume of primitive cell of 
the metal layer. The consequence of this effective field is to induce, in linear response, a 
deviation in longitudinal magnetization 

SM Z = Xzz SH eff (B3) 

where Xzz is the longitudinal susceptibility for interface spins of the FI layer and SM Z is 
related to the magnon density 5M z (x) = —2fi B 5n m (x). Thus, the relation between spin 
accumulation of non-equilibrium electron density and nonequilibrium magnon density at 
the metal-FI interface is 

- = £5n m (B4) 

where 

9 & (B5) 



JsdD (s F )xzzVo' 

with M.. A/To — 

(2tt) 



Xzz can be derived in linear response by Xzz = ^ijt\ H ^ with M z = Ms — / dqN^ (q). 
In the limit Dq^^ 3> /c#T, we find 

= VjM^L Lil _ (e-A./.BT) (B6) 
4 (irDf 2 2 V ' 

where Li s (z) = YlkLi F ^ s the polylogarithm. Inserting Eq.(B6) into (B5), we get the final 
result e of Eq. (18). 
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